Density of states in coupled chains with off-diagonal disorder 
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We compute the density of states p(e) in A'' coupled chains with random hopping. At zero energy, 
p(e) shows a singularity that strongly depends on the parity of A'^. For odd N, p{e) oc l/|eln''£|, 
with and without time-reversal symmetry. For even N , p{e) oc | Ine] in the presence of time-reversal 
symmetry, while there is a pseudogap, p{e) oc |£lne|, in the absence of time-reversal symmetry. 
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Since the works of Dyson, it has been known that one 
dimensional systems with off-diagonal disorder may show 
anomalous behavior in their density of states and their 
localization properties 0. This anomalous behavior has 
reappeared in many reincarnations, ranging from Dyson's 
original application of masses coupled with springs, to the 
random hopping model |^,|| , quantum XY spin chains , 
supersymmetric quantum mechanics Q , one-dimensional 
classical diffusion in a random medium Q, narrow-gap 
semiconductors 0, and charge density wave materials 
In all these systems, it has been shown that both 
the density of states and the localization length diverge 
at zero energy. These two divergencies are closely linked 
by a theorem by Thouless |^ . 

Recently, it was found that the localization length of a 
wire consisting of A^ of such chains shows a remarkable 
dependence on the parity of A^ ]To| , pd| : At zero energy, i.e. 
at the center of the band, the localization length diverges 
for odd A^ (odd A^ includes the pure one-dimensional case 
A^ = 1), while wavefunctions remain localized for even A^. 
Moreover, it was found that staggering in the hopping 
parameter (as the consequence of a Peierls instability), 
may lead to an additional set of delocalized states for 
all A^ > 1, regardless of parity jl^. For comparison, 
in the pure one-dimensional case A^ = 1 staggering al- 
ways enhances localization. In this paper we consider 
the density of states (d.o.s.) for the A^-chain wire with 
off-diagonal disorder. Unlike in the one-dimensional case, 
this is a problem that requires separate attention, since 
the Thouless theorem which links d.o.s. and localization 
length does not hold for A^ > 1 coupled chains. More- 
over this problem is also of relevance for the random flux 
model, which is a special example of off-diagonal ran- 
domness, and whose localization properties and density 
of states near e = are the subject of an ongoing debate 
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To be specific, we consider the d.o.s. for the 
Schrodinger equation 
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where '0„ is an A^-component wavefunction and t„ an 
NxN hopping matrix. The index n = 1, . . . , L labels the 



site index along the wire. The length of the wire, mea- 
sured in units of the lattice spacing, is L. The reason why 
this system can display anomalous behavior near zero en- 
ergy is the existence of a particle-hole or chiral symmetry, 
which is absent in the presence of on-site disorder: Under 
a mapping (— 1)"'0„, the Hamiltonian TC changes 

sign, Ti — > —H. As a result, the eigenvalues of H occur 
in pairs ±£. There are two mechanisms by which the 
chiral symmetry is known to affect the d.o.s. near e = 0. 
First, level repulsion of the eigenvalue e with its mir- 
ror image — e causes a universal suppression of the d.o.s. 
near e = This suppression appears on the scale 

of a level spacing, is independent of the geometry, but 
becomes unimportant in the thermodynamic limit. The 
second mechanism, for which several different descrip- 
tions exist 1^-^,01 , is special for a (quasi) one-dimensional 
geometry, and survives in the thermodynamic limit. It 
is responsible for the divergence of the d.o.s. in the pure 
one-dimensional case. Below we present a calculation of 
the d.o.s. in the multi-chain case, combining ideas from 
the Fokker-Planck approach to localization in multichan- 
nel quantum wires and the calculation of the d.o.s. 
in the one-dimensional random hopping model |^,^ . 

We first state our main results. We find that the parity 
dependence that was previously obtained for the localiza- 
tion properties is also present in the density of states. For 
odd A^, the d.o.s. diverges at zero energy according to 
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where u; is a dimensionless parameter governing the ran- 
domness in the t„. The bandwidth is chosen as the unit 
of energy. In Eq. ^ and in the remainder of the paper 
we assume e > 0. The d.o.s. for £ < then follows from 
p{—e) — p{£). The form of the divergence (||) is indepen- 
dent of whether time- reversal symmetry is broken or not, 
i.e. of whether the hopping matrices t„ are generically 
real or complex. For even N , in contrast, the density 
of states strongly depends on the presence or absence of 
time-reversal symmetry (labeled by the parameter (3=1 
or 2, respectively), 



p{£) CX |e/w;2|^~^|ln(wV£)|, A^ even. 
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In the presence of time-reversal symmetry, pie) shows a 
logarithmic divergence as e — > 0, while in the absence 
of time-reversal symmetry a pseudo-gap is opened, p{e) 
vanishes at e = 0. Such a strong dependence on time- 
reversal symmetry is a remarkable result. An effect of 
comparable magnitude appears in the suppression of the 
gap by a weak magnetic field in a normal metal in the 
proximity of a superconductor. In the remainder of this 
paper we derive the results (||) and (^) assuming a specific 
statistical model for the hopping matrices i„ in Eq. (|^). 
The details of this model are not relevant, the singular 
behavior of the d.o.s. near zero energy being governed by 
the fundamental symmetries of the Hamiltonian Ti only. 
As an illustration of the general validity of our result, we 
close with a comparison to numerical simulations. 

As we are interested in the d.o.s. of the random hop- 
ping model (|^) in the thermodynamic limit, the bound- 
ary conditions at the two ends of the chain are not im- 
portant. For convenience we choose hard wall boundary 
conditions, — = 0, L being the length of the chain. 
We can solve Eq. recursively in terms of a sequence 
of hermitian N x N matrices a„, 

a„Vn+i = tii'ri, an = "t^ (s + a„_i)"^ (4) 

The boundary condition at n = implies ao — 0. Evalu- 
ating the Schrodinger equation at n — L then yields that 
e is an eigenvalue of Ti. if and only if 

eV'L = -t^_i'(/'L-l = -UL-l'ipL, (5) 

i.e. if ol-i has an eigenvalue — e, or alternatively, in view 
of Eq. (^), if fli has a diverging eigenvalue. 

As a statistical model for the i„ that contains all the 
relevant symmetries, we parameterize the t„ in terms of 
the generators 

t„ = e^", (6) 

where Wn is a real (complex) matrix for P — I (2), and 
choose the matrices Wn from independent Gaussian dis- 
tributions with mean and variance given by 

{{Wn)^AWn);,) = ^w'i3[s^,6,^ - '-g^s,jp^], 

{(Wn),.) = i(-l)"A5^.. (7) 

Here r/ governs the fluctuations of tr Wn fl^ and A mea- 
sures the staggering of the hopping parameter. (A is the 
gap size that the staggering would induce in the absence 
of disorder.) We assume that the disorder and staggering 
are weak, and that the energy is small compared to the 
bandwidth (w^. A, e < 1). 

The matrix has eigenvalues a^, that can be param- 
eterized as = tan((^/2) (/i = 1, . . . ,N). As we have 
discussed below Eq. (0), the energy e is an eigenvalue 
of the Hamiltonian, if and only if there is an angle (/>^ 
with 0^ = TT. For general e, however, none of the 0^ will 



be equal to tt. Nevertheless, we can use the angles (/)^ to 
compute the (disorder averaged) density of states. Hereto 
we first note that aL+2 = in the absence of disorder, 
staggering, and for e = 0. Then, taking disorder, stag- 
gering, and a finite energy into account, and considering 
the length L as a fictitious "time" , the angles perform 
a Brownian motion on the unit circle, which is such that 
upon increasing L, they move around the circle in posi- 
tive direction. The rate at which the i/i^ pass through tt 
as we increase L, i.e. their current, equals the number of 
states per unit length N(s) with energy between and 
£. (This is a generalization of the node-counting theorem 
|l6| used to compute the d.o.s. for = 1 [||j^.) For com- 
parison, we remark that, in the absence of disorder, the 
angles (/)^ move around at a constant speed oc e, resulting 
in a constant d.o.s. With disorder, their motion acquires 
a random (Brownian) component, which dramatically af- 
fects their average speed, and hence the density of states, 
as we shall see below. 

For a quantitative description a different parameteriza- 
tion of the eigenvalues proves to be more convenient, 

a^=tan(0^/2) = e"-. (8) 

The variables are restricted to the two branches 
Iniu^ = and Imu^ = tt in the complex plane, see 
Fig. We refer to these as lower and upper branches, 
respectively. Noting that the are related to the angles 
(j)fj, on the unit circle, we see that a (fictitious) particle 
with coordinate that vanishes on one of the branches 
at ±oo reappears at the opposite branch, as indicated by 
the arrows in Fig. Upon increasing L by two, the 
change according to ^ + (5w^, where, to lowest 
order in w, e, and A, the average and variance of the 
increments 5u^ are 

(Suf,) = 2g cosh + 2A + ^ coth ~ "'^ ^ 

{Su^Su,) = Aw^ [5^, - (1 - ri)/N] . (9) 

Taking L as a continuous variable, their distribution 
function P{ui, . . . ,u^; L) obeys the Fokker-Planck equa- 
tion 




J =Y[smh^[{u^-u,)/2]. (10) 

For large L the solution of Eq. (|l^) acquires a steady 
state carrying a current j(s) equal to the integrated d.o.s. 
N(e) in the thermodynamic limit. Unfortunately, except 



2 



t 1 

max " I '^max 

FIG. 1. (a) Two branches for the coordinate u. (b) 
In the simplified model, the branches are truncated at 
i^^max = ln(ui^/e). 

in the case iV = 1, where a solution in closed form is pos- 
sible 01 , it is notoriously problematic to find the steady- 
state solution of a Fokker-Planck equation of the type 
(|l0|), due to the lack of detailed balance Instead, 
below we present a qualitative analysis of the Brownian 
motion process described by Eq. (^0|), illustrating the 
mechanisms that lead to the anomalous behavior of the 
d.o.s. near e = 0. We give a detailed account of the case 
rj — 1, followed by a brief discussion of the general case 
?7 7^ 1 and a comparison with numerical simulations. 

Let us first identify the relevant parameters in the 
Brownian motion process ([l0|). There are N fictitious 
Brownian "particles" with coordinates {fJ- = 1, ■ • ■ , N) 
and diffusion coefficient D = w^. Three forces F^, Fa, 
and act on the particles, arising from the presence of 
disorder, staggering, and energy, respectively. The force 
Fw = coth[(u^ — Ui,)/2] is a repulsive interaction, 

with a hard-core on the same branch and a soft core on 
different branches. (We adopt a convention where the 
particles have unit mobility.) Staggering causes a con- 
stant force field Fa = A that favors motion to the right 
on both branches. Finally, F^ = ztecoshu pushes the 
particles to the left (right) on the lower (upper) branch 
and thus causes the nonzero steady- state current. For 
small energies e <C w^, motion is governed by the force 
Fg for large |Rew| ^ Kmax, where 

Umax = ln(w^/e), (11) 

while diffusion and the forces F^ and F^ are dominant 
for small u, \Keu\ ^ Umax- Following Rcf. we now 
approximate our model by truncating the branches at 
|Reu| = Umax, see Fig. |]b, and adding a one-way move 
towards the upper (lower) branches at the end-points, see 
Fig. |]b. These simplifications don't alter the functional 
dependencies of p{e) for e — > 0, though they may affect 
prefactors. 

For = 1, the "particle" with coordinate u needs a 
time (2ui„ax)^/2£' to diffuse to the end point of a branch, 
as one can verify from a solution of the one-dimensional 
diffusion equation on the line —Umax < u < Umax, with 
hard wall boundary conditions at u = —Umax and absorb- 
ing wall boundary conditions at m = Umax |^ • Hence the 
current is j{e) = N{e) oc D/Au^^^, and after differenti- 
ation with respect to e, one finds Eq. (0) for the density 
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FIG. 2. Schematic picture of the forces on the fictitious 
particles in the simplified model of Fig. For N = 2, their 
mutual repulsion slows the particles down and traps them 
near the end-points (a). When they eventually meet, the re- 
pulsion enhances travel (b), and quickly restores the situation 
(a). The case A*' = 3 is similar to A*' = 1: Two particles are 
trapped near the ends as a result of their repulsion, while 
there is no net force on the third particle (c). For A'' = 4, 
all particles are trapped near the end points (d). The repul- 
sive force on the two middle particles U2 and U3 is smallest, 
and these two dominate the current according to a scenario 
similar to the case N = 2 (see a and b). 

of states. Diffusion effectively speeds up the particle, 
explaining the enhancement of the d.o.s. relative to the 
clean case ||ll-§,0. For N = 2, the picture is completely 
different. As a result of their mutual repulsion, the two 
particles with coordinates ui and U2 get trapped near the 
end points, say at ui — —Umax and U2 = Umax + iT^- Now 
the particles have to diffuse out of their traps against 
their repulsive interaction, until they eventually meet at 
Re Ml = Re U2 and the repulsive force F^ starts to favor 
travel (see Figs. ||a and b). Such a process costs a large 
time, which can be calculated from the diffusion equation 
for two particles on a line with hard wall boundary condi- 
tions at Re Ml = —Umax and Reu2 = Umax and absorbing 
boundary conditions at Reui = Reu2. We find a cur- 
rent j{e) = (4umaxF3/D2)exp(-2F„Mmax/-D), resulting 
in the d.o.s. (^). (The prefactor Umax arises from the de- 
generacy of the meeting point on the line.) We conclude 
that the d.o.s. for two coupled chains has only a loga- 
rithmic singularity for real hopping disorder, and a pseu- 
dogap for complex disorder. The strong /3-dependence of 
the d.o.s. stems from the /3-dependence of the interaction 
force F^,. 

The qualitative behavior of the d.o.s. for general N > 2 
depends crucially on the parity of A^ and closely resem- 
bles the scenarios we have outlined above for A" = 1 and 
N — 2. If Af is even, all particles get "trapped" near the 
ends of the branches, half of them on the lower branch 
near —Umax, and half of them on the upper branch near 
Mmax, like in the case A^ = 2, cf. Fig. The repulsive 
interaction force F.^ is smallest for the two particles that 
are closest to the origin. These two particles dominate 
the current, resulting in a d.o.s. of the form (^). If A^ is 
odd, on the other hand, the picture is like that of the case 
A^ = 1, see Fig. ||c. All particles get trapped at the two 
ends, except for one "free" particle, for which diffusion is 
not slowed down by the interaction forces (the repulsive 



forces from the other particles cancel exactly). Hence, 
for odd N, the d.o.s. is of the form (||). This even-odd 
effect is reminiscent of that found for the conductance at 
e = 0, where the existence of a delocalized state could be 
attributed to the existence of a similar "free" transmis- 
sion eigenvalue . 

The effect of staggering A is to add a constant force 
pointed to the right on both branches. As a function of 
the staggering strength, the system alternates between 
behavior corresponding to even and odd N. Repeat- 
ing the above analysis, we find that with staggering, the 
d.o.s. shows the maximum (|^) if A = {N + 1 — 2j)it;^/3/2, 
J = l,...,iV, whereas it shows the minimum if 
A = (N - 2j)w'^p/2, j = 1, . . . , TV - 1. For all other 
values of A, for e — > we have 
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(12) 



The case of arbitrary 77 7^ 1 is not much different from 
the case 77 = 1 we considered above; up to prefactors 
the e-dependence of the d.o.s. is not changed. The only 
exception is the case 77 = 0, = 2, when the center 
of mass ui + U2 is pinned. As a result, the degeneracy 
giving rise to the logarithm in Eq. (||) is lifted, and the 
logarithmic prefactor vanishes. 

We conclude with a comparison to numerical simula- 
tions for the d.o.s. in a quantum wire on a square lattice 
with a width between A^ = 1 and A = 4, and a length 
L between 10"^ and 10^. For (3 — 1 and also for A^ = 1 
the hopping amplitudes are taken from a uniform dis- 
tribution in the interval [0.5,1.5], while for (3 = 2 the 
random flux model |12| is used, where the randomness 
is introduced only via the random phases of the hopping 
amplitudes. Results of an average over 4 x 10^-10^ disor- 
der realizations are shown in Fig. ^ The agreement with 
our theoretical results, Eqs. (||) and (||), is excellent. 
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for A = 1). 
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